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ABSTRACT 

Narrow-band imaging surveys aimed at detecting the faint emission from the 5007 
A [O III] line of intracluster planetary nebulae in Virgo also probe high redshift z ~ 3.1 
Lya emitters. Here we report on the spectroscopic identification of 9 Lya emitters at 
z = 3.13 with fluxes between 2 x 10~ 17 and 2 x 10~ 16 erg cm -2 s -1 obtained with the 
FORS spectrograph at Unit 1 of the ESO Very Large Telescope (VLT UT1). 
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The spectra of these high redshift objects show a narrow, isolated Lya emission with 
very faint (frequently undetected) continuum, indicating a large equivalent width. No 
other features are visible in our spectra. Our Lya emitters are quite similar to those 
found by Hu (1998), Cowie & Hu (1998) and Hu et al. (1998). For a flat universe with 
H =70 km s _1 Mpc -1 and q =0.5 (1^=0), the Lya luminosity of the brightest source 
is 1.7 x 10 9 L0, and the comoving space density of the Lya emitters in the searched 
volume is 5 x 1CT 3 Mpc~ 3 . 

Using simple population synthesis models, on the assumption that these sources are 
regions of star formation, we conclude that the nebulae are nearly optically thick and 
must have a very low dust content, in order to explain the high observed Lya equivalent 
widths. For the cosmological and star formation parameters we adopted, the total stellar 
mass produced would seem to correspond to the formation of rather small galaxies, some 
of which are perhaps destined to merge. However, one of our sources might become a 
serious candidate for a proto-giant spheroidal galaxy if we assumed continuous star 
formation, a low mass cutoff of 0.1 M© in the IMF, and a flat accelerating universe 
with ^o=0.2 and ft A =0.8. 

The implied star formation density in our sampled comoving volume is probably 
somewhat smaller than, but of the same order of magnitude as the star formation 
density at z ~ 3 derived by other authors from Lyman-break galaxy surveys. This 
result agrees with the expectation that the Lya emitters are a low-metallicity (or low- 
dust) tail in a distribution of star forming regions at high redshifts. Finally, the Lya 
emitters may contribute as many H-ionizing photons as QSOs at z ~ 3. They are 
therefore potentially significant for the ionization budget of the early universe. 

Subject headings: cosmology: observations — early universe — galaxies: clusters: in- 
dividual (Virgo cluster) — galaxies: intergalactic medium — galaxies: evolution - 
galaxies: formation 



1. Introduction 

Lya emitting galaxies at high redshift are of much interest for studies of galaxy formation. 
While early surveys failed to find such objects (e.g. Thompson et al. 1995), recent narrow-band 
imaging searches with spectroscopic follow-up identified a number of Lya emitters at high z, first in 
fields near high-redshift QSOs (Hu & McMahon 1996) and later in blank fields (Cowie & Hu 1998; 
Hu et al. 1998; Pascarelle et al. 1998). These Lya emitters have typically very faint continua and 



1 Based on observations collected at the European Southern Observatory, Cerro Paranal, Chile; ESO programmes 
63.N-0530 and 63.1-0007 
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high Lya equivalent widths, i.e., may represent an early phase of galaxy formation when substantial 
amounts of dust had not yet formed. 

Here we report on the discovery of 9 Lya emitters at redshift z = 3.1 during our program of 
studying the intracluster planetary nebula (PN) population in the Virgo cluster, which uses similar 
narrow-band imaging techniques to identify candidate PNs. The first indication of the existence 
of a diffuse intracluster stellar population in Virgo was the discovery by Arnaboldi et al. (1996) 
that a few PNs in the galaxy NGC 4406 (M 86; redshift -230 km s^ 1 ) have redshifts typical of the 
Virgo cluster (around +1300 km s _1 ). Subsequently, direct evidence for red giant stars belonging 
to this stellar population was reported by Ferguson et al. (1998). Because PNs offer the chance to 
measure radial velocities and perhaps even abundances for such a diffuse population, a search for 
intracluster PNs in different positions across the Virgo cluster started immediately and produced 
several dozens of PN candidates (Mendez et al. 1997, Feldmeier et al. 1998) in a total surveyed 
area of 0.23 deg 2 . 

The "on-band/off-band" narrow-band filter technique used to discover the PN candidates (see 
e.g. Jacoby et al. 1992) allows the detection of a single emission line. This is not necessarily the 
desired [O III] A5007; it might be another emission line at higher z, redshifted into the on-band 
filter. For example [O II] A3727 at z = 0.35, or Lya at z = 3.13. In previous work (Mendez et al. 
1997) we argued that most of our detections had to be real PNs because (a) the surface density of 
emission line galaxies derived from earlier studies was not high enough to explain all the detections; 
(b) the luminosity function of the detected sources is in good agreement with the PN luminosity 
functions derived in several Virgo galaxies (Jacoby et al. 1990). 

We have used the first ESO Very Large Telescope (VLT) unit (UT1) with the Focal Reducer and 
Spectrograph (FORS) in multi-object spectroscopic mode, and the 4-m Anglo- Australian Telescope 
(AAT) with the 2-degree-field (2df) fiber spectrograph. Our purpose was to confirm the nature of 
the Virgo intracluster PN candidates, to measure their radial velocities, and (in the VLT case) to 
detect the faint diagnostic lines required for abundance determinations. 

The present paper reports on the results of the VLT+FORS observations. The AAT+2df 
observations, which confirm the existence of Virgo intracluster PNs through the detection of both 
the A4959 and A5007 [O III] lines (Freeman et al. 1999), will be presented and discussed by Freeman 
et al. (2000, in preparation). 

2. Observations and first results 

For the purpose of detecting faint nebular diagnostic lines with VLT+FORS (ESO programme 
63.N-0530), we selected Field 1 of Feldmeier et al. (1998) as our first priority, because it had 
the brightest PN candidates, with a luminosity function cutoff about half a magnitude brighter 
("^5007 = 25.8) than elsewhere in the Virgo cluster (see Feldmeier et al. 1998). The total area 
of Field 1 is 256 arcmin 2 . We also wanted to use FORS (ESO programme 63.1-0007) to verify 
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the PN nature of candidates in the smaller "La Palma Field" (50 arcmin 2 , Mendez et al. 1997), 
with magnitudes 777,5007 between 26.8 and 28.6. The lack of bright PNs in the La Palma Field was 
understandable in Mendez et al. (1997) as a consequence of the sample size effect: if the total PN 
sample is small, the chance of finding a bright PN is too small. 

The observations were made with FORS at the VLT UT1 on the nights of 11/12, 15/16, 16/17 
and 19/20 April 1999. Since the FORS field is slightly smaller than 7x7 arc mins (with the standard 
collimator), we selected a portion of Field 1 where several of the brightest candidates were located, 
and took on-band and off-band images of the selected portion of Field 1, and of the La Palma Field, 
on the night 11/12, using FORS in imaging mode. 

The on-band and off-band interference filters we used have, respectively, central wavelengths 
5039 and 5300 A, and FWHMs 52 and 250 A. Some 50 stellar images in the short on-band and off- 
band FORS exposures (10 min and 3 min, respectively) were used together with the corresponding 
stellar images in the Kitt Peak and La Palma discovery images to define coordinate transformations 
that gave the pixel values of the positions of the PN candidates in the FORS images, given the 
pixel values of their positions in the discovery images. In this way it was possible to define the 
FORS slitlet positions with sufficient accuracy (in fact, the brightest PN candidates were visible 
in the short on-band FORS exposures, allowing us to verify directly the accuracy of the pixel 
transformation, with typical errors below 0.5 pixels, which is equivalent to 0.1 arc sec). 

On the night of 15/16 April, R.P.K. and R.H.M. started FORS multi-object spectroscopy 
(MOS) in Field 1, using grism 300V without any order separation filter, for maximum spectral 
coverage. This grism gives a spectral resolution of about 10 A, at 5000 A, for slitlets 1 arc sec 
wide. Of the 19 available slitlets, 10 were placed at the positions of PN candidates (numbers 1, 4, 
5, 6, 16, 17, 21, 26, 27, 42, ordered by brightness as measured by Feldmeier). These candidates had 
magnitudes 777.5007 between 25.8 and 26.4. The remaining 9 slitlets were placed on stars or galaxies 
in the field, in order to check the slitlet positioning (which was done by taking a short exposure 
without grism through the slitlets) and to help locate the dispersion lines as a function of position 
across the field, which is important in the case of spectra consisting of isolated emission lines. 

After taking 5 exposures of 40 min each in Field 1, all with the same slitlet configuration, and 
having made preliminary on-line reductions, it was clear that the [O III] emission line candidates 
were not PNs. Object 1 showed one strong isolated emission line at 5021 A. It cannot be [O III] 
A5007 because there is no hint of the companion line A4959 at the corresponding wavelength. Object 
5 was identified as a starburst with z = 0.35; it shows narrow emission lines with a continuum, 
A3727 is redshifted into the on-band filter, and FL3 and the [O III] lines 4959 and 5007 are visible, 
also redshifted with z = 0.35. Object 17 was not detected. The remaining 7 objects were identified 
as continuum objects. Many of them are galaxies: they show redshifted emissions. 

In summary, none of the candidates tested in Field 1 with VLT+FORS are PNs. This result 
in fact solves a problem, because these Field 1 candidates were surprisingly abundant and were 
somewhat brighter than typical PNs in the Virgo galaxies (see the discussion by Feldmeier et 
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al. 1998). The high percentage of continuum objects (7 out of 9 detections) indicates that the 
off-band exposure by Feldmeier et al. was not deep enough. A re-examination of the Field 1 
images by J. Feldmeier has subsequently confirmed that due to sudden changes in the seeing and 
transparency the off-band does go to a brighter limiting magnitude than intended. Briefly, the Field 
1 exposures consisted of three 3600s on-band exposures and five 600s off-band exposures taken at 
the Mayall 4-m telescope. The additional off-band exposures were intended to compensate for the 
transparency, which was decreasing at the time of the exposures. Unfortunately, although these 
additional exposures do partially compensate for the change in transparency, the seeing increased 
by 0.3 arcsec, and consequently the mean off-band exposure is not deep enough. 

On the other hand, Freeman et al. (1999) did confirm spectroscopically some fainter PN 
candidates in Field 1. Taken together, these results imply that the surface density of the intracluster 
stellar population originally estimated for Field 1 has to be reduced (Freeman et al. 2000, in 
preparation). 

In view of the result in Field 1, we immediately pointed the VLT UT1 to the La Palma field. 
In this field 12 PN candidates had been found, 11 reported in Mendez et al. (1997) and one found 
afterwards. The distribution of the 12 PN candidates in the sky made it impossible to assign slitlets 
to all of them in one slitlet configuration. In our first configuration we defined 1 arc sec wide slitlets 
for 9 PN candidates and 2 objects suspected to be QSOs or starbursts because they were visible 
(although much weaker) in the off-band discovery image (Mendez et al. 1997). On the April nights 
15/16 and 16/17 we completed 5 MOS exposures (40 min each) of the initial slitlet configuration, 
and on April 19/20 we took 3 additional MOS exposures of the La Palma Field (again 40 min each) 
with a different slitlet configuration, which allowed to add one PN candidate not observed before. 

Thus spectra for a total of 10 PN candidates were acquired in the La Palma Field, and 7 were 
detected. They all show an isolated and narrow emission at wavelengths from 5007 to 5042 A. In 
all cases this is the only feature visible in the spectrum. The other 3 PN candidates, which are 
relatively faint, were not detected. Perhaps their slitlets were slightly misplaced, although we would 
think the errors in the pixel transformation were too small to have any effects on the detectability. 
Perhaps some of these sources, if they are not PNs, have variable brightness. 

Of the 2 QSO or starburst candidates, one was confirmed as a QSO at z=3.13, showing a 
typical broad-lined Lya and C IV A1550. The other one appears to be a starburst because it shows 
one strong, isolated and narrow emission line, with faint continuum. 

The success rate for emission line detection in the La Palma Field was satisfactorily high. It 
may be useful to give a few numbers for comparison with other searches, complementing information 
given in Mendez et al. (1997). The La Palma on-band image had a limiting magnitude msoo7 = 
28.7, equivalent to a flux of 10~ 17 erg cm" 2 s _1 . We would collect the same amount of photons 
through the on-band filter from a star of visual magnitude 25.5. The off-band image had a limiting 
magnitude 0.2 mag fainter. The search for candidates was done by blinking the on-band versus the 
off-band image. To compare with the selection criteria used e.g. by Steidel et al. (1999b) for their 



-6- 



narrow-band imaging survey, we define on-band and ofF-band magnitudes so that they coincide for 
an average star. All our emission-line candidates, being fainter or absent in the off-band image, 
have positive colors offband— onband. We can give only lower limits to the colors of candidates 
undetected in the off-band image. All tested emission-line candidates with colors above 1.0 mag 
have been spectroscopically confirmed. The lower limits to the colors of the 3 objects that were 
not detected with FORS are below 1.0 mag. Our 2 QSO or starburst candidates, detected in both 
images, have colors of 2.1 and 1.7 mag, respectively. 

3. Analysis of the La Palma Field spectra 

The CCD reductions were made using IRAF 4 standard tasks. After bias subtraction, flat field 
correction and image combination to eliminate cosmic ray events, the object spectra were extracted 
and the sky background subtracted. Then the He-Ar-Hg comparison spectra were extracted and 
the object spectra were wavelength calibrated. Spectrograms of the standard stars G138-31 and 
G24-9 (Oke 1990) were used for the flux calibration. 

We have designated the La Palma Field sources with LPF plus a number ordered according 
to the brightness in the discovery image. LPFnew is the latest PN candidate, found after the 
discovery paper (Mendez et al. 1997) was published. Object LPFsl is the QSO, and LPFs2 is the 
object identified as a starburst from the start, because of its stronger continuum. Fl-1 and Fl-5 
are objects 1 and 5 in Feldmeier's Field 1. Fl-1 turns out to also have a very faint continuum, 
barely visible in the final processed spectrogram. 

Figs. 1 to 3 show, respectively, the spectra of the QSO, the starburst with visible continuum, 
and one of the La Palma Field PN candidates showing no detectable continuum. All the LPF PN 
candidates look very similar, with only one emission line detected across the whole spectrum. 

Having found no direct evidence of [O III] A4959, which would be expected if the detected 
emission line were A5007, we can put an upper limit to the percentage of our emission-line objects 
that can be PNs. After rejecting objects that show a continuum, which clearly cannot be PNs, we 
proceed in the following way: 

(1) shift the spectra, so that the wavelengths of the detected emission lines fall at 5007 A. 

(2) normalize the intensities of the emission lines to the same value, e.g. 300. 

(3) add all the spectra and measure the intensity of the resulting A4959. If it is 50, for example, 
comparing to the expected value of 100 (since A5007 was defined to be 300) we can argue that 50% 
of the objects must be PNs. 



4 IRAF is distributed by the National Optical Astronomical Observatories, operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under contract to the National Science Foundation of the U.S.A. 
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The result of this test is shown in Fig. 4, where we have added the normalized spectra of the 
7 PN candidates detected in the La Palma Field. The complete absence of A4959 indicates that, at 
most, one of the 7 candidates can be a PN. This conclusion is based on the noise level, not on any 
marginal detection of A4959. Besides, there is in Fig. 4 a hint of a weak continuum, not detectable 
in the individual spectra, which reinforces the rejection of these candidates as PNs. 

4. Identification of the detected emission line as Lya 

Having rejected [O III]A5007, because A4959 is not visible in Fig. 4, we consider the alterna- 
tives: 

(1) [O II] A3727 at z = 0.35 was confirmed in one case (object 5 in Field 1) but can be rejected 
in all other cases because we do not see H/3, [O III] AA4959, 5007 and Ha at the corresponding 
redshifted wavelengths. This is illustrated in Figs. 5 and 6. 

(2) Mg II A2798 at z = 0.79 can also be rejected for a similar reason: in this case we do not see 
[O II] A3727 at the expected redshifted wavelength, as illustrated in Fig. 7. The same argument 
can be applied to other lines: assuming C III A1909, we do not see A2798; and so on. 

We conclude that the isolated emission line must be Lya at z = 3.1. This identification is 
supported by the strength of the line: since we see at most a very faint continuum, the equivalent 
width is fairly large, typically 200 A (observed) and 50 A (rest frame). We have mentioned in 
Mendez et al. (1997) that very few starburst galaxies show, for example, [O II] A3727 stronger 
than 100 A in equivalent width. 

5. Implications for the surface density of intracluster PNs in Virgo 

A reliable estimate of the surface density of intracluster PNs in Virgo will have to await a 
survey of a sufficiently large area on the sky, which is currently in progress. Our results to date 
and some preliminary conclusions are the following. 

No PN candidates have been confirmed in the La Palma Field. Only 5 of the 12 candidates 
have a chance of remaining as PNs: 2 were not tested and 3 were tested but not detected. Since 
we have identified most of these candidates as Lya emitters, it is clear that the intracluster PN 
sample size in the La Palma Field and the inferred surface density must be substantially smaller 
than we estimated. 

On the other hand, the AAT+2df multi-object fiber spectroscopy has confirmed the existence 
of intracluster PNs in Fields 1 and 3 of Feldmeier et al. (1998). This will be reported in detail by 
Freeman et al. (2000, in preparation). These PNs are brighter than rri5oo7 = 27. Freeman et al. 
will show that the contamination by Lya emitters at these brighter magnitudes is not as important 
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as in the La Palma Field. Thus it appears that the fraction of Lya emitters in [0 III] narrow-band 
selected samples is magnitude-dependent, increasing towards fainter values of 7715007. 

The lack of bright PNs in the La Palma Field implies a lower surface density than in other 
Virgo cluster positions, and may indicate some degree of clustering in the distribution of the diffuse 
intracluster population. 

Our VLT+FORS observations have shown that a spectroscopic confirmation of intracluster PN 
candidates, involving the detection of both the A4959 and A5007 [O III] emission lines, is necessary. 

6. Observed properties of the high-redshift Lya emitters 

In Table 1 we have collected some basic information about the narrow-lined sources: the 
measured Lya fluxes, measured Lya wavelengths, redshifts, and Lya equivalent widths (Wa in 
almost all cases lower limits, because the continuum is not detected). The measured wavelengths 
provide yet another argument favoring the interpretation of all these sources as unrelated to the 
Virgo cluster: they are randomly distributed across the on-band filter transmission curve, with no 
concentration at the Virgo cluster redshift (between 5020 and 5030 A). See Fig. 8. 

Notice that object LPF3 is brighter than measured in the discovery image. The reason is that 
in this case Lya falls near the edge of the on-band filter transmission curve. LPF3 is remarkable 
also for the very large lower limit to its Lya equivalent width. 

The observed fluxes are between 2 x 10~ 17 and 2 x 10~ 16 erg cm~ 2 s _1 . We are looking into a 
small redshift range, defined by the transmission of the on-band filter used at La Palma, of about 
Az = 0.04. The La Palma Field covers 50 arcmin 2 . The total Lya flux measured within our 
"discovery box" is 5 x 10~ 16 erg cm~ 2 s _1 . Adopting z = 3.13, Ho=70 km s _1 Mpc -1 and qo=0.5 
we get a luminosity distance of 1.8 xlO 4 Mpc, which implies Lya luminosities for our sources 
between 2 x 10 8 and 2 x 10 9 L©. The total Lya luminosity within the sampled volume is 5 x 10 9 
LQ. The sampled comoving volume is 1650 Mpc 3 , which gives from 8 sources a comoving space 
density 5 x 10~ 3 Mpc -3 . 

These numbers depend on our assumptions about the cosmological parameters: for a flat 
universe with f2o=0.2, f2A = 0.8, and the same z and Ho as above, the luminosity distance becomes 
3 x 10 4 Mpc, the sampled comoving volume 10 4 Mpc 3 , and the Lya luminosities become larger by 
a factor 2.8. 

7. Starburst models and derived quantities 

What produces the narrow Lya emission? Given rest frame equivalent widths below 200 A, 
the most probable source of ionization is massive star formation (Chariot and Fall 1993). Active 
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Table 1. Observed and derived properties of Lya emitters. 



object 


flux a 


A(A) 


z 


w(A) b 


M stars c 


SFR d 


M s &SFR e 


LPFs2 


17 


5011 


3.121 


33 


21-420 


7-140 


5-10 


LPF1 


7 


5011 


3.121 


>55 


9-75 


3-25 


2-3 


LPF2 


6 


5026 


3.133 


>32 


7-150 


3-50 


2-4 


LPF3 


10 


5042 


3.146 


>175 


8-11 


3-4 




LPF4 


3 


5007 


3.118 


>14 


4-270 


1-90 


1-5 


LPF6 


2 


5035 


3.140 


>70 


3-14 


1-5 


0.4-0.5 


LPF8 


2 


5034 


3.140 


>37 


3-44 


1-15 


0.6-1 


LPFnew 


2 


5010 


3.120 


>25 


3-65 


1-22 


0.6-2 


Fl-1 


15 


5021 


3.129 


45 


19-200 


6-67 


4-6 



a measured Lya fluxes in units of 10 -17 erg cm~ 2 s^ 1 . Using the lu- 
minosity distance of 1.8 xlO 4 Mpc, this column also gives the Lya lu- 
minosities, in units of 10 8 L©. 

b equivalent widths of Lya transformed into the rest frame (the mea- 
sured widths are z + 1 times larger) . 

c Possible range of total stellar mass formed, in units of 10 6 M©, as- 
suming starbursts of 3 x 10 6 yr duration; see section 7 and Fig. 14. 

d Possible range of star formation rates, in M© yr -1 , assuming star- 
bursts; see section 7. 

e Possible range of total stellar mass formed, in units of 10 9 M©, as- 
suming continuous star formation over 10 9 yr; see section 7 and Fig. 15. 
The same numbers also represent the corresponding star formation rates 
in M© yr -1 . 
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galactic nuclei (AGNs) might be an alternative, although there is no hint of C IV 1550 in our 
spectra, as shown in Fig. 9. We have taken massive star formation as our working hypothesis. 

We have made a simple population synthesis model to explore some basic properties of the 
stellar population responsible for the ionization of the H II regions. Our main interest in doing 
this analysis is to estimate star formation rates and densities; we would like to verify if these Lya 
sources correspond to the formation of massive giant galaxies or rather to the formation of smaller 
structures. 

We adopt a standard initial mass function (IMF), f(M) oc M -2 ' 35 , stellar evolutionary models 
for low metallicity (Z=0.001; this choice of metallicity will be justified in next paragraph) from 
Schaller et al. (1992), and ionizing fluxes from NLTE model atmospheres with winds (Pauldrach 
et al. 1998), again for a low metallicity, in this case that of the SMC (5 times below solar). 

The relation between the number of stellar Lyman continuum photons Ni y c and the Lya 
luminosity can be obtained from a simple recombination model: 

L(Lya) = hz/(Lya) 0.68 X B N LyC (1) 

where 0.68 is the fraction of recombinations that yield Ly a (Case B, see e.g. Storey and Hummer 
1995), and Xb is the product of the fraction of Nl v c really absorbed times the fraction of Lya 
photons that really escape. We necessarily have < Xb < 1- Note that Xb ~ 1 requires both an 
optically thick nebula and very low dust content, because the large number of Lya scatterings in 
an optically thick nebula will lead to their absorption by dust grains, if such grains are present in 
any significant number. This explains our choice of a low metallicity. How low is low? Since we do 
not know much about dust properties, we have used empirical information provided by Chariot and 
Fall (1993) in their Fig. 8, which shows Lya equivalent widths as a function of oxygen abundance 
in nearby star-forming galaxies. In that figure we find that Lya equivalent widths larger than 
20 and 50 A are associated respectively with oxygen abundances below 25% and 10% solar. In 
future work, to obtain more quantitative constraints, we intend to carry out Lya radiative transfer 
calculations in the presence of dust and velocity fields as an extension of the work by Hummer and 
Kunasz (1980), Hummer and Storey (1992) and Neufeld (1990, 1991). 

Figs. 10 and 11 show Lya equivalent widths for single stars as a function of Xb and of the 
stellar T c g. We need both a large Xb and many massive, hot main sequence stars in order to 
produce the observed Lya equivalent widths. 

We have explored two different histories of star formation in a low- metallicity stellar population: 
(a) a starburst, and (b) continuous star formation. These two alternatives are defined as star 
formation extending over a time (a) similar to (3 x 10 6 years) and (b) much longer than (10 9 years) 
the duration of a main sequence OB star. 

Figs. 12 and 13 show the resulting contour plots of Lya equivalent widths as a function of 
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Xb and of the maximum main sequence mass in the population. In these figures, the quantities 
on the x axis allow different interpretations. In the case of a starburst, smaller values of M max can 
be interpreted to represent an increasing age of the starburst, or in other words the time elapsed 
since the starburst happened; as the starburst grows older, the most massive stars are removed 
from the main sequence. In the case of continuous star formation, the value of M max indicates at 
which mass the integration of the IMF was stopped. In other words, the left part of the plot shows 
a case in which very massive main sequence stars were not formed. 

We can obtain Lya equivalent widths > 50 A only for Xb > 0.5 (continuous star formation) 
or Xb > 0.3 (starburst). This points to almost completely optically thick, extremely dust-poor 
nebulae. It is easy to understand why a large Xb is necessary. If it is small, many stellar ionizing 
photons are lost. In order to explain the observed Lya fluxes we must add more massive stars, but 
these stars make a strong contribution to the continuum, and therefore the Lya equivalent width 
must decrease. 

It may seem surprising that we can get large Lya equivalent widths at so low values of M max . 
The reason is that the Schaller et al. main sequence at low metallicity is shifted to rather high T c g 
because of the low stellar opacity. This means that lower mass objects on the main sequence give 
much more ionizing flux than at solar metallicity. 

The source LPF3 can be explained only as a very young starburst with Xb ~ 1, because of its 
very large Lya equivalent width. 

We set the lower mass limit of the Salpeter IMF at 0.5 MQ. For a maximum stellar mass of 120 
MO, this gives an average mass of 1.65 MQ, which is then the conversion factor between number 
of stars and total stellar mass. The average mass decreases only slightly if we decrease M max , and 
only if -Mmax is interpreted as in Fig. 13, because in that case very massive stars are not formed. 
The average mass does not decrease if M max is low due to the age of the starburst, because the 
very massive stars are assumed to have formed and evolved away from the main sequence. 

Figs. 14 and 15 show contour plots of the logarithms of Lya luminosities in a plane where 
the x axis is the same M max used in Figs. 12 and 13, and the y axis is the product of Xb times 
the total number of stars produced over the total duration of the star formation (3 x 10 6 and 10 9 
yr, respectively). The number of stars required to produce a given Lya luminosity depends on the 
value of M max . Since we cannot determine M max empirically or derive it from first principles, we 
have considered the full range of M max , from 120 MQ to as low as permitted by the limits on the 
Lya equivalent widths; sometimes as low as 20 MQ. This produces rather large uncertainties in 
the derived masses, particularly in the case of starbursts. A comparison of Figs. 14 and 15 also 
shows that, as expected, in the case of continuous star formation many more stars are needed to 
obtain a given Lya luminosity. 

The results of these mass estimates are listed in Table 1. The total masses of stars and star 
formation rates turn out to be very sensitive to the star formation history. For continuous star 
formation and young starbursts, i.e. high values of M max , our strongest sources have SFRs of the 
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order of 10 M© yr . However these numbers increase dramatically if we consider older starbursts, 
and we cannot rule out SFRs higher than 200 M© yr -1 in a few cases. 

Adding up the masses and SFRs in Table 1, and using the average stellar mass of 1.65 M©, 
we have built Table 2, which provides the total number and total mass of stars formed, and the 
resulting star formation rate, needed to explain our total Lya luminosity of 5 x 10 9 L© in the La 
Palma Field, for the cases of starbursts and continuous star formation. Here we have not added 
the numbers for Fl-1, which belongs to another field. 

How uncertain are these numbers? First of all, the real SFR might be higher than we inferred 
for our sources because of extinction by dust. However, we consider this to be unlikely because it 
would be necessary to argue that there is dust in the foreground (to produce extinction) but not 
inside (it would destroy the Lya emission). We do not argue in terms of low metallicity because 
there is at least one case of a very metal-deficient blue compact dwarf galaxy, SBS 0335-052, where 
dust patches are clearly visible (Thuan and Izotov 1999). Note however that this source has Lya in 
absorption, implying that extinction is accompanied by destruction of Lya emission, as expected. 
Based on this reasoning we do not include a correction for extinction in our Tables. 

The reader may argue that geometry probably plays the dominant role in dictating how many 
Lya photons escape, making our attempt to estimate star formation rates a futile exercise. However, 
our sources are a rather special case. Geometry is useful when we need to explain why a star-forming 
region shows Lya in absorption; see e.g. Kunth et al. (1998). But our sources do show strong Lya 
emission, and furthermore, show a very faint continuum. The absence of continuum (i.e. the large 
equivalent width of Lya) is very important because it precludes the existence of large amounts 
of hot stars. If there is any significant destruction of Lya photons, then, in order to explain the 
observed Lya fluxes we must add more massive stars. But these additional stars contribute to the 
continuum, and the Lya equivalent width decreases too much. Thus the absence of continuum 
acts as a "safety valve" restricting higher values of the star formation rate. As explained above, 
extinction by dust is not very likely. 

On the other hand, the derived masses and SFRs are strongly dependent on the lower mass 
limit of the IMF. If we selected a lower mass limit of 0.1 instead of 0.5 M© there would be no 
change in the number of stars needed to produce the ionizing photons, but we would be adding a 
lot of low-mass stars; the total mass estimate for a given total Lya luminosity would increase by 
a factor 1.9. The SFRs would have to be corrected by the same factor. The ranges of masses and 
SFRs in Tables 1 and 2 do not include this source of uncertainty, but we will take it into account 
in the discussion. 

The adopted cosmological parameters also have some influence. We have adopted Ho=70 km 
s _1 Mpc -1 and qo=0.5, which implies a sampled comoving volume of 1650 Mpc 3 . If we take, for 
example, a total SFR of 50 M© yr -1 , we get a star formation density in the comoving volume of 
0.03 M© yr -1 Mpc -3 . If we adopted a universe with Qq=0.2, Q^=0.8, but with the same z and 
Hq, then our star masses and SFRs would be larger by the same factor 2.8 as for the luminosities. 
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However the star formation density mentioned above would drop from 0.03 to 0.014 M© yr 
Mpc~ 3 because of the enlarged sampled comoving volume of 10 4 Mpc 3 . 

8. Discussion 

Recent blank field searches for Lya emitters have been successful at redshifts from 2.4 to 5: 
Hu 1998, Cowie and Hu 1998, Pascarelle et al. 1998, Hu et al. 1998. The sources that have been 
verified spectroscopically (Hu 1998, Hu et al. 1998) are like ours: strong, narrow, isolated emission 
line identified as Lya, with a similar range of equivalent widths. 

One may ask why the earlier searches failed to detect such emission-line sources. For example 
Thompson et al. (1995) reported no detection in a total area of 180 arcmin 2 down to a l-a limiting 
flux 10 times fainter than the flux of our brightest source, LPFs2. This might be attributed to 
a very low surface or space density in their directions, but since they looked in many directions 
this interpretation looks improbable. The surface density in the La Palma Field does not seem 
to be abnormally high. The 8 sources we detected in the 50 arcmin 2 of the La Palma Field and 
within our redshift range of 0.04 are equivalent to 14400 emitters deg -2 per unit z with Lya fluxes 
above 1.5 x 10~ 17 erg cm~ 2 s -1 . This is a lower limit to the true density of such sources; if some 
of the three photometric candidates not detected and of the two candidates not tested are Lya 
emitters with similar fluxes, the true density could be higher by up to 60%. The density inferred 
from the present spectroscopic sample is similar to the 15000 deg~ 2 per unit z reported by Hu et 
al. (1998) from Keck searches in different regions of the sky, with the same limiting flux. Our 
sources are brighter than those detected with the Keck telescope (Hu 1998, Cowie and Hu 1998, 
Hu et al. 1998), but our redshift is smaller; see Fig. 16. In summary, the surface density of Lya 
emitters in the La Palma Field is of the same order of magnitude as for the sources detected by Hu 
et al. Pascarelle et al. (1998) report an order-of-magnitude range of space densities at z = 2.4, in 
some cases even higher than ours. More surveys will probably clarify whether or not there is any 
structure or clustering in the distribution of Lya emitters. 

Concerning masses and star formation rates, if we adopt continuous star formation or young 
starbursts, where star formation has not yet (or has just) stopped, then from Table 1 we read that 
our strongest sources show a SFR not higher than about 10 M© yr -1 , similar to the maximum SFRs 
reported by Hu et al. (1998), and short of what would be expected from a proto-giant spheroidal 
galaxy forming more than 10 11 M© in 10 9 years. We cannot completely rule out the possible 
existence of an inconspicuous stellar population, already formed, which would be detectable only 
in the infrared, but we would find it difficult to explain the very low dust content in that case. The 
possible existence of undetected neutral and molecular H gas would provide additional means to 
scatter and eventually destroy the Lya photons (e.g. Neufeld 1990), and therefore we consider it 
unlikely. Thus from the available evidence we would seem to be witnessing the formation of small 
subgalaxies, some of which are perhaps destined to merge. This conclusion is also supported by 
the total mass produced, even in the case of continuous star formation, and would be unchanged if 
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we adopted the lowest IMF mass limit of 0.1 M©. If we also adopted the flat accelerating universe 
with £lo=0.2 and fiA=0.8, which implies larger distances, luminosities and SFRs, the starbursts 
would still point to small entities, but the continuous star formation in the case of LPFs2 would 
be able to produce several times 10 10 M0. 

On the other hand, continuous star formation implies the previous production of supernovae 
and metals, which makes it more difficult to explain a very low dust content. For that reason we 
consider young starbursts more likely than continuous star formation or old starbursts. But in 
dealing with this subject we prefer to be cautious and leave all conceivable options open. 

What would happen if we allowed for older starbursts, in which star formation has stopped 
10 7 years ago? Then the mass of stars formed would increase substantially, but notice that in such 
cases the mass produced after star formation has stopped is not larger than 10 9 M©. Therefore 
even in this extreme case starbursts are not able to make a proto-giant spheroidal galaxy out of 
any of our sources. This conclusion would not be weakened by assuming the accelerating universe 
mentioned above. 

Now we try to estimate lower and upper limits for the star formation density in the sampled 
comoving volume that corresponds to the La Palma Field, assuming that we can rule out production 
of Lya photons through AGN activity. A lower limit for the SFR of 15 MQ yr _1 comes from the 
continuous star formation case in Table 2. The maximum SFR in Table 2 (old starbursts) is rather 
improbable, because it is obtained assuming that all the starbursts ended at the same time, some 
10 7 years before the Lya photons we detected were emitted. A more reasonable upper limit is 
obtained assuming that one of the starbursts is at the right age to produce the maximum SFR, 
while all the others make smaller contributions. Let us assume that the source produced by the 
only allowed old starburst is LPFs2, which would then contribute 140 M© yr" 1 . Adding the almost 
negligible contributions from the other sources, we obtain in this way a range of plausible total 
SFRs in the comoving volume between 15 and 170 MQ yr _1 . Allowing now for the uncertainty in 
the lower mass limit of the IMF, we get a range between 15 and 320 MQ yr~ x , which implies a 
star formation density between 0.009 and 0.19 MQ yr _1 Mpc -3 . For comparison, Hu et al. (1998) 
obtained, from their own samples, 0.01 in the same units. 

In order to compare our star formation density with other available information, we adopt in 
this paragraph Ho=50 km s _1 Mpc -1 and qo=0.5, which is what other authors have done, e.g. 
Steidel et al. (1999a). This gives luminosities twice as large as in our earlier choice, and a sampled 
comoving volume of 4500 Mpc 3 . Our resultant star formation density drops slightly to between 
0.007 and 0.14 MQ yr -1 Mpc~ 3 . This is plotted in Fig. 17 together with other data collected from 
several galaxy surveys by Steidel et al. (1999a). Our star formation density in Lya emitters is 
comparable to that in other star forming sources at that redshift. Our sources are probably the 
low-metallicity (or low-dust) tail in a distribution of star forming regions at high redshifts. This is 
underlined by the fact that in our sampled volume we also have a QSO (LPFsl) with apparently 
higher metal abundances, judging from the strength of the C, N, O lines in its spectrum. 
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As already remarked by Hu et al. (1998), since we expect lower metallicity at higher redshifts, 
we should expect the strong Lya emitters to become more frequent at higher redshifts relative to 
Lyman-break galaxies, which normally have weak or absent Lya emission and are therefore presum- 
ably more metal-rich. Note, however, that low metallicity does not necessarily imply emission in 
Lya: the blue compact dwarf galaxy SBS 0335-052 (Thuan and Izotov 1999), with Z = Z Q /41 and 
Lya in absorption, provides a beautiful cautionary note. Kunth et al. (1998) have argued that the 
geometry and velocity structure of the interstellar medium play an important role in determining 
the strength of the Lya emission, and this may complicate the comparison between Lya emitters 
and Lyman-break galaxies. 

Finally, we consider our Lya emitters as possible sources of ionization of the intergalactic 
medium at high redshift (see e.g. Madau et al. 1998). Keeping the same cosmological parameters as 
in last paragraph, we get a total L(Lya), in our sampled comoving volume, of 10 10 L©. Converting 
this into Nl v c using Eq. (1), on the assumption that Xb = 0.5, we get 7 x 10 54 photons s -1 . 
Since the sampled comoving volume is 4500 Mpc 3 , this implies a ionizing photon density of 1.5 
xlO 51 photons s _1 Mpc~ 3 . This is the total number of ionizing photons produced by the stars. It 
is very difficult to estimate what fraction is available for ionization of the intergalactic medium, 
because it depends on the relative contributions of the two factors that enter into Xb- In at least 
one case (LPF3) we have Xb ~ 1, which means that this source cannot contribute many ionizing 
photons. Assuming optimistically that there remain 10 51 ionizing photons s _1 Mpc~ 3 available, 
our Lya emitters seem to contribute about 0.5 times as much as the star-forming galaxies at z = 3 
found in Lyman-break galaxy surveys (see e.g. Fig. 2 of Madau et al. 1998). We conclude that the 
contribution by Lya emitters may be comparable in order of magnitude to what QSOs provide at 
z = 3. 

9. Summary of conclusions and perspectives 

We have discovered a population of high-redshift Lya emitters in our sample of Virgo intr- 
acluster PN candidates obtained with an "on-band/off-band" filter technique. Our VLT+FORS 
spectra show that the Lya emitters at z = 3.13 look very similar to those discovered in other fields 
at other redshifts by Hu et al. (1998). Only a narrow and strong Lya emission is visible in their 
spectra; no other spectral line, and no (or at most a very weak) continuum. On the assumption 
that Lya emission is produced by massive star formation, we have estimated the total mass of 
stars formed and star formation rates, and we have estimated the star formation density in our 
sampled comoving volume. The Lya emitting nebulae must be nearly optically thick and extremely 
dust-poor, probably indicating a very low metallicity. The total mass formed and the SFRs appear 
to suggest that we are witnessing the formation of rather small galaxies. This conclusion depends 
on several assumptions about the IMF, star formation history and cosmological parameters. There 
is one source (LPFs2) that might qualify as a proto-giant spheroidal if we assumed continuous star 
formation, a lower mass cutoff of 0.1 M© in the IMF, and a flat accelerating universe with Qq=0.2 
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and J7a=0.8. However, to assume continuous star formation implies the previous production of 
supernovae and metals, making it more difficult to explain the very low dust content required by 
the observed spectra of our sources. We are more probably dealing with young starbursts. 

Taking all sources into account, the implied star formation density in our sampled comoving 
volume is probably somewhat smaller than, but of the same order of magnitude as the star formation 
density at z ~ 3 derived by other authors from Lyman-break galaxy surveys. This result agrees 
with the expectation that the Lya emitters are a low-metallicity (or low-dust) tail in a distribution 
of star forming regions at high redshifts. Finally, the Lya emitters may contribute as many H- 
ionizing photons as QSOs at z ~ 3. They are therefore potentially significant for the ionization 
budget of the early universe. 

More extensive surveys at different redshifts will be needed to build a luminosity function 
for the Lya emitters and to decide if they show any evidence of clustering and evolution as a 
function of redshift. HST images might be able to resolve these sources; their morphology might 
offer clues about their star formation processes. High-resolution spectroscopy of Lya would provide 
important kinematic information about the interstellar medium in these galaxies. We need infrared 
spectra to try to detect forbidden lines and either determine or put some upper limits to the 
metallicity. It might also be possible to clarify to what extent the production of Lya photons can 
be attributed to AGN activity. The firm detection of an infrared continuum would help to constrain 
the characteristics and total mass of the stellar populations through comparison with population 
synthesis models. 

The implications of our results for the surface density of the diffuse intracluster stellar popula- 
tion in the Virgo cluster are not clear yet. Future searches for intracluster PNs will have to include 
the spectroscopic confirmation of the candidates, by detection of the two bright [O III] emissions 
at 4959 and 5007 A. From the spectroscopic work of our group to date, including the confirmation 
of 23 intracluster PNs by the detection of the two [OIII] lines by Freeman et al. (1999), the fraction 
of high-redshift Lya sources in the on-band/off-band samples appears to be higher at faint magni- 
tudes. In the specific case of the "La Palma Field" we have not found any intracluster PNs, perhaps 
suggesting some degree of dumpiness in the distribution of the intracluster stellar population. In 
other Virgo fields with brighter PN candidates, the fraction of high-redshift Lya emitters among 
the detected sources appears to be about 25% (see Freeman et al. 2000, in preparation, for a more 
extensive discussion). 

RPK would like to thank the director and staff of Steward Observatory, Tucson, for their 
hospitality and support during an inspiring sabbatical which among other results has led to the 
ESO proposal 63.N-0530. Thanks also to Stefan Wagner for helpful comments. RPK and RHM 
express their gratitude to the ESO staff at Cerro Paranal, who contributed with their efforts to 
make this VLT run an enjoyable experience. 
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Fig. 1. — The spectrum of QSO candidate LPFsl. The strong, broad emissions (Lya + N V A1240 
and C IV A1550) confirm that this is a QSO at a redshift z = 3.13 (the observed wavelength of 
Lya is 5024 A). A faint continuum is visible only longward of 3760 A, which is consistent with the 
position of the Lyman break. Other features visible are O VI A1035 at 4275 A, O IV + Si IV A1402 
at 5790 A, and He II A1640 at 6775 A. 

Fig. 2. — The isolated emission in the spectrum of starburst candidate LPFs2. A continuum is 
clearly seen, as expected from the off-band detection in the discovery paper. 

Fig. 3. — The isolated emission in the spectrum of PN candidate LPF3. Here no continuum is 
detected, consistent with the off-band non-detection in the discovery paper. 

Fig. 4. — Addition of the normalized spectra of 7 PN candidates (no individual continua detected). 
The normalization procedure is described in section 3. The arrow indicates where [O III] A4959 
should be, if the visible emission were [O III] A5007. The absence of 4959 leads to conclude that 
probably none, and at most one, of the 7 candidates can be a PN. Note also that in this added 
spectrum a very faint continuum may be visible. 

Fig. 5. — The upper spectrum corresponds to an anonymous galaxy in the La Palma Field which 
happened to show strong emission at [O II] A3727, H/3 and [O III] AA4959, 5007. We use this 
object as reference; its continuum has been arbitrarily rectified for easier comparison. Below this 
star-forming galaxy we show the spectrum of Fl-5 (object 5 in Field 1), redshift-corrected so that 
the strong emission line detected with the on-band filter falls at 3727 A. Fl-5 shows the same set of 
starburst emission lines as the reference galaxy, and therefore its redshift is confirmed to be 0.35. 
Below we have plotted the spectra of our sources, all redshift-corrected in the same way as for 
Fl-5. The levels of zero intensity are separated by 500 counts. None of these sources shows other 
emission lines at the relevant wavelengths. 

Fig. 6. — The upper spectrum corresponds to another anonymous galaxy, this time in Field 1, 
which shows emission at A3727, A5007 and Ha. The arrow indicates the position of Ha. Below 
it, from top to bottom, we show the spectra of 6 of our sources, all redshift-corrected as described 
in Fig.5. The levels of zero intensity are separated by 500 counts. None of these sources shows 
any convincing evidence of emission at the expected wavelength of Ha. The "inverted P Cyg 
features" in the spectra of LPFs2 and LPF1 are artifacts produced by imperfect sky subtraction 
in the presence of very strong sky emission lines. We have verified by careful inspection of the 
combined CCD spectrograms that no real features are visible at those wavelengths. The spectra of 
the sources which are not shown here do not reach so large wavelengths, because of their positions in 
their respective fields. This figure and the previous one show that we cannot attribute the detected 
emission lines to [O II] A3727. 

Fig. 7. — Here all spectra have been redshift-corrected so that the strong emission line detected 
with the on-band filter falls at Mg II 2798 A. The levels of zero intensity are separated by 500 
counts. Nothing is visible at the expected position of [O II] A3727, indicated with an arrow. 
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Fig. 8. — The full line is the transmission curve of the on-band filter used at La Palma. The 
histogram shows the wavelength distribution of the emission lines. One object corresponds to 1 
unit along the y axis. Our sources are more or less uniformly distributed across the filter curve. 

Fig. 9. — Spectra of our sources in the rest frame of Lya. At the top is the QSO LPFsl. The 
levels of zero intensity are separated by 500 counts. Only the QSO shows C IV A1550. At some 
wavelengths a slightly larger noise level is noticeable; this is due to imperfect sky subtraction in 
the presence of very strong sky emission lines. In several cases where the disturbance by the sky 
lines was too bad, we have replaced the spectrum by a straight line at zero intensity. 

Fig. 10. — Equivalent width of Lya emission produced by a single star, as a function of the stellar 
T e ff, calculated for Xb = 0.3. 

Fig. 11. — Equivalent width of Lya emission produced by a single star, in the plane T e g vs. Xb- 
The contour lines are labeled with the Lya equivalent widths in A. 

Fig. 12. — Contour plots of Lya equivalent widths in A, as a function of Xb and M max (the 
maximum main sequence mass in the population), for the case of a starburst. See text. 

Fig. 13. — Same as Fig. 12, for the case of continuous star formation. See text. 

Fig. 14. — Contour plots of the logarithms of Lya luminosities, in units of Lq, for the case of a 
starburst. Given a Lya luminosity, it is possible to read, for different values of M max , the product 
of Xb times the total number of stars. 

Fig. 15. — Same as Fig. 14, for the case of continuous star formation. 

Fig. 16. — The logarithms of observed Lya fluxes (in erg cm~ 2 s -1 ) as a function of redshift. The 
plus signs are our sources, and the vertical bars represent the ranges of fluxes observed by Hu 
(1998) at redshifts 3.4 and 4.6. The dashed lines represent the observed fluxes for Lya luminosities 
of 10 42 and 10 43 erg s _1 , if we assume Ho=70 km s _1 Mpc -1 and qo=0.5. 

Fig. 17. — Logarithms of star formation density in M© yr -1 Mpc -3 , as a function of redshift. The 
lines show the results derived from galaxy surveys, collected by Steidel et al. (1999a): neglecting 
(dashed) and correcting for (solid) extinction by dust. These curves are derived assuming continuous 
star formation. Our star formation densities at z = 3.13 have been adapted, for consistency, 
to the same cosmological parameters used by Steidel et al. The square and diamond represent, 
respectively, our results assuming continuous star formation and starbursts. Both symbols have 
been slightly displaced horizontally for legibility. 
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Table 2. Star formation needed to explain the total Lya luminosity in the La Palma Field. 
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Fig. 2 and 3. — Figures 2 and 3 




Fig. 4. — Figure 4 
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Fig. 5. — Figure 5 
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Fig. 12 and 13. — Figures 12 (starburst) and 13 (continuous star formation) 




Fig. 14 and 15. — Figures 14 (starburst) and 15 (continuous star formation) 
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Figures 16 and 17 



